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I .  Introduction 

Experimental  data  on  relaxation  phenomena  in  diverse  areas  of  condensed 
matter  physics  are  quite  generally  found  to  exhibit  slower  than  exponential 
decay  for  long  times  in  the  form  (see  ref.  1  and  references  therein) 

exp[-atb],  a>0,  0<b<l.  (1) 

The  physical  origin  of  this  behavior  will  be  discussed  here  from  several 
points  of  view.  These  include  the  necessity  of  time  dependent  transition 
rates  (TDTR)  in  relaxation  theory  and  their  treatment  via  a  time  scale  trans¬ 
formation,  the  constraints  on  the  form  of  the  TDTR  from  quantum  mechanical 
considerations  and  a  theorem  in  Fourier  transforms,  the  requirements  on  micro¬ 
scopic  distributions  in  frequency  of  the  underlying  relaxing  entities  such  as 
dipoles  or  charge  carriers  from  basic  probability  theory  of  limit  distribu¬ 
tions  for  sums  of  independent  random  variables,  and  consideration  of  the 
generalized  master  equation  (GME)  with  an  underlying  nonMarkoffian  stochastic 
process  with  TDTR,  all  leading  to  long  time  nonexponential  relaxation  given  by 
Eq.  (1).  We  also  discuss  the  relationship  of  Eq.  (1)  with  the  mathematical 
framework  of  the  continuous  time  random  walk  (CTRW),  GME,  and  the  H- theorem  of 
statistical  mechanics.  A  quantum  mechanical  model  is  described  which  leads  to 
Eq..(l)  and  in  addition  predicts  a  relaxation  time  renormalization  which 
alters  activation  energies,  and  temperature  and  molecular  weight  dependences. 
These  predictions  have  been  verified  experimentally  in  electronic  materials, 
ionic  conductors,  glasses,  polymer  melts,  and  other  materials.  The  calcula- 
tional  ease  of  the  time  scale  transformation  method  is  demonstrated  for  an 
important  physical  example.  The  general  applicability  of  Eq.  (1)  to  elucidate 
the  vast  array  of  experimental  results  in  relaxation  phenomenon  in  condensed 
matter  is  briefly  described.  Thus  TDTR  leading  to  Eq.  (1)  is  found  to  be 
essential  in  describing  relaxation  phenomena. 
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II.  Time  Scale  Transformation 


The  traditional  phenomenological  description  of  relaxation  is  based  on 
stochastic  processes  with  time  independent  transition  rates  (TITR)  leading  to 
the  usual  exponential  decay.  In  view  of  the  experimentally  observed  nonexpo¬ 
nential  behavior,  given,  by  Eq.  (1),  it  is  necessary  to  consider  TDTR  in  the 

2  t>“l 

stochastic  theory.  We  have  shown  that  a  special  class  of  TDTR  given  byjt  , 

can  be  related  to  TITR  by  means  of  a  time  scale  transformation,  for  which  the 
relaxation  function  is  of  the  form  Eq.  (1)  and  the  time  scale  transformation 
is  the  monomial,  8(t)  =  at**.  This  comes  ab^ut  by  the  requirement  that  the 
TDTR  be  such  that  the  relaxation  function  associated  with  it  depend  only  on 
the  ratio  t/T  where  t  is  a  constant  characteristic  time  and  in  the  time  trans¬ 
formed  frame  this  process  becomes  one  with  TITR  with  corresponding  relaxation 

function  exp(-0/T  )  where  T  is  another  constant.  The  only  time  transforma- 
s  s 

tion  that  has  the  property  0(t)  =  which  converts  a  TDTR  to  a  TIDR  as  above 
is  the  monomial. 

III.  The  Paley-Wiener  Theorem  and  TDTR 

The  transition  rate  Can  be  calculated  quantum  mechanically  and  the  use  of 
the  golden  rule  leads  to  TITR.  However  the  golden  rule  is  an  approximation 
and  this  breaks  down  whenever  the  energy  difference  of  the  two  states  between 
which  transitions  are  taking  place  goes  to  zero.  In  the  relaxation  regime, 
typically  in  the  range  of  101®  Hz  and  below,  such  characteristic  energy  dif¬ 
ferences  may  approach  zero  and  a  more  careful  calculation  leads  to  TDTR.  A 
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general  bound  on  TDTR  is  provided  by  the  Paley-Wiener  theorem  in  Fourier 
transform  theory,  which  is  found  to  lead  to  a  lower  bound  on  the  relaxation 
function  in  the  form  Eq.  (1).  The  requirements  of  the  Paley.-Wiener  theorem  in 
our  context  are  semiboundedness  of  the  energy  spectrum  of  the  Hamiltonian  and 
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the  quadratic  integrability  of  the  relaxing  states,  both  of  which  are  physi¬ 
cally  acceptable. 

The  relaxing  system  is  a  complex  many-body  system  whose  low  energy  exci¬ 
tations  below  approximately  lO1^  Hz  determine  TDTR.  A  detailed  knowledge  of 
the  many-body  excitations  is  not  available  at  present.  However,  the  many-body 
density  of  states  of  a  random  matrix  Hamiltonian  has  been  deduced  by  Wigner 
sometime  ago  and  this  can  be  used  here  to  obtain  TDTR.  Interestingly,  this 
also  leads  to  Eq.  (1).  More  details  of  this  dynamical  model  are  given  in  Sec. 
VIII. 

IV.  Limit  Distributions  of  Low  Energy  Excitations 

Given  the  universal  form  of  the  macroscopic  relaxation  function  in 
Eq.  (1),  one  is  led  to  consider  a  microscopic  mechanism  underlying  relaxation 
phenomenon  in  general.  We  may  think  of  the  relaxing  body  as  being  composed  of 
some  species  of  relaxing  entities  such  as  electric  dipoles  in  a  dielectric, 
charge  carriers  in  a  semiconductor,  etc.,  each  of  which  has  identically  dis¬ 
tributed  energy  variables  with  the  same  energy  distribution.  The  macroscopic 
energy  distribution  function  p(e)  associated  with  the  relaxing  body  for  low 
energies  is  the  limit  distribution  of  normalized  sums  of  the  microscopic 
energy  variables.  A  similar  procedure  is  involved  in  statistical  mechanics 
for  deriving  equilibrium  ensemble  distributions.  The  resulting  macroscopic 
density  distribution  p(s)  necessarily  obeys  the  relationship^ 

p(ae+b)  =  Cajajj/aJj^dc'pCajCc-c' )+b)p(a26'+b2)  (2) 

for  every  set  of  parameters  ai  >  0,  b^  with  the  corresponding  a  >  0,  b  related 
to  the  parameters  a^,  b^.  The  relaxation  function  is  proportional  to  the 
modulus  square  of  the  Fourier  transform,  | c(t)J  * ,  of  the  energy  distribution 
function  p(e).  For  distributions  obeying  Eq.  (2)  this  is  necessarily  of  the 
form^ 
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Jc(t)J  2  =  ,  a>0,  0<b<2  (3) 
If  in  addition,  the  spectrtim  of  the  Hamiltonian  is  semibounded  as  in  Sec.  Ill, 
then  the  Paley-Wiener  theorem  applies  and  the  values  of  b  have  to  be  in  the 
smaller  range  0<b<l,  leading  us  again  to  Eq.  (1).  Foi.  such  microscopic 
models,  we  see  that  the  Paley-Wiener  bound  becomes  exact. 

V.  An  Oscillator  Relaxation  Model 

A  chain  of  coupled  oscillators  can  serve  as  a  model  of  a  Brownian 
particle  in  contact  with  a  heat  bath.**  For  arbitrary  oscillator  interactions, 
the  Brownian  particle  obeys  a  generalized  Langevin  equation  with  a  time 
dependent  friction  coefficient,  and  a  Gaussian  random  force  determined  solely 
by  the  Hamiltonian  of  the  coupled  oscillators  and  the  canonical  distribution 
of  the  initial  coordinates  and  momenta.  The  Paley-Wiener  theorem  can  be 
applied  to  this  problem  if  the  normalized  frequency  spectrum  of  the  oscillator 
assembly  is  bounded  from  below.  The  Fourier  transform  of  the  spectrum  is  the 
momentum  autocorrelation  function  of  the  single  Brownian  oscillator  and  this 
is  then  bounded  by  a  function  of  the  form  Eq.  (1).  If  the  bound  is  assumed,  a 
time  scale  transformation  can  be  used  to  convert  the  generalized  Langevin 
equation  into  the  classical  Langevin  equation  with  a  time- independent  friction 
coefficient  and  a  Gaussian  white  noise  random  force.  Transforming  back  to  the 
t-frame,  we  find  that  the  covariance  of  the  external  random  force  takes  the 
form  of  white  noise  modulated  by  a  monomial  friction  coefficient,  t^  *.  It 
can  be  noted  that  a  covariance  of  the  same  form  is  obtained  in  the  long  time 
limit  from  fractional  Brownian  motion. ^  However,  this  connection  to  frac¬ 
tional  Brownian  motion  is  rather  tenuous  because  the  form  of  the  noise  does 
not  uniquely  specify  the  stochastic  equation  of  motion. 
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VI.  Continuous  Tine  Random  Walk  and  Generalized  Master  Equation  * 

There  exists  a  mathematical  framework  of  the  continuous  time  random  walk 

8  9 

(CTRW)  and  the  associated  generalized  master  equation  (GME)  for  studying 

time-dependent  problems.  The  basic  quantity  entering  the  formalism  is  the 

waiting  time  distribution  function  «Ji(t) ,  which  is  the  probability  density 

function  for  the  time  between  the  arrival  of  a  walker  at  a  given  lattice  point 

and  the  initiation  of  the  next  step  to  another  site.  If  the  waiting  time 

distribution  function  is  taken  as  minus  the  derivative  of  Eq.  (1),  then  the 

CTRW  framework  can  be  applied  to  relaxation  problems  in  condensed  matter. 

When  the  moments  u  of  the  waiting  time  distribution  functions  are  all  finite, 
m 

the  Markoffian  master  equation  is  an  appropriate  description  for  times  which 

are  large  compared  with  t  =sup(u  /m!)1'  .  For  «|i(t)  given  by  minus  the  deriva- 

m 

* 

tive  of  Eq.  (1)  it  is  found  that  t  is  finite  for  l<b<2  and  infinite  for 
0<b<l.10  Thus  a  non-Markoff ian  description  is  always  required  for  relaxation 
functions  of  the  form  Eq.  (1). 

In  most  applications  of  CTRW,  in  order  to  have  a  <{»(t)  which  is  universal 
for  all  sites,  one  has  to  -perform  an  average  over  all  possible  configurations 
of  the  walkers.  This  physical  picture  is  inappropriate  for  materials  where 
short  range  order  exists,  which  is  the  common  situation  in  condensed  matter. 
A  model,  described  in  Sec.  VIII,  in  which  a  relaxation  function  of  the  form 
Eq.  (1)  arises  from  a  fundamental  mechanism  that  is  operative  at  every  hop  of 
the  walker  from  any  lattice  site,  produces  a  «j>(t)  universal  for  all  sites 
without  averaging. 

VII.  CTRW,  GME,  and  the  H-Theorem 

In  CTRW,  the  weighting  time  distribution  function  is  related  to  the 

q 

memory  function  of  the  GME.  If  the  GME  is  to  describe  a  time- dependent 
statistical  mechanical  phenomenon,  it  must  obey  the  H- theorem,  which  is  a 
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statement  of  the  second  law  of  thermodynamics.  The  H-theorem  cannot  accommo¬ 
date  arbitrary  memory  functions  and  therefore  also  arbitrary  waiting  time 
distribution  functions.  The  choice  of  minus  the  derivative  of  Eq.  (1)  for  the 
waiting  time  distribution  in  the  CTRW  along  with  the  time  scale  transformation 

of  Sec.  II  enables  us  to  demonstrate  the  validity  of  the  H-theorem  for  this 

.  .  10 
choice. 


VIII.  Model  for  Low-Energy  Excitations 

A  model*  in  which  a  relaxing  primary  species,  such  as  a  reorienting 
electric  dipole,  a  hopping  charge  carrier  or  a  flexing  polymer  chain,  inter¬ 
acts  with  a  new  class  of  very  low  energy  states  (correlated  states  or  CS) 
produces  a  relaxation  function  of  the  form  Eq.  (1)  as  well  as  correlates  the 
dispersion  with  shifts  in  activation  energies,  molecular  weight  dependence, 
temperature  dependence,  and  other  quantities.  A  random  matrix  formulation  of 
energy  level  distributions  can  be  used  to  study  the  properties  of  the  CS.  It 
is  found  that  for  materials  having  short-range  spatial  order,  the  effect  of 

such  CS  excitations  is  to  alter  a  constant  TITR  x  *  to  a  TDTR  given  by  t  * 

o  o 

exp(-ny)  (u>ct)  n  for  times  greater  than  approximately  u>c*.  Here  n  is  a 
measure  of  the  coupling  strength  of  the  CS  and  is  restricted  to  0<n<l.  The 
cutoff  frequency  u>c  depends  on  material  structure  but  is  typically  of  the 
order  10*°  Hz  for  glasses,  polymers,  and  amorphous  metals  and  semiconductors, 


and  yS0.577. 

The  above  TDTR  results  in  a  relaxation  function  of  the  form  Eq.  (1)  with 
b=l-n  and  a  shifted  relaxation  time  Tp  s  a  given  by 

t,  =  I(l-»)enl'  <*> 


The  CS  model  is  therefore  able  to  identify  the  physical  meaning  of  the  param¬ 
eters  in  Eq.  (1)  and  is  able  to  make  predictions  other  than  dispersion. 


The  implications  of  Eq.  (A)  are  widespread.  For  example  if  t  is  temper¬ 
ature  activated  with  activation  energy  E^,  then  the  relaxation  process  will 
appear  to  be  activated  with  energy  E*=E^/(l-n).  This  prediction  has  been 
verified  in  detail  for  transient  current  in  chalcogenide  glasses. **  The 
correct  value  of  E*  is  predicted  from  the  measured  value  of  E^  and  the  value 
of  n  obtained  from  the  dispersion  of  the  current.  In  addition,  the  values  of 
Ejt  correlate  correctly  with  changes  in  the  dispersion  as  the  value  of  n  is 
changed  by  doping. 

Another  result  of  Eq.  (4)  is  the  prediction  of  molecular  weight  depen- 

12 

dences  of  viscoelastic  quantities  for  polymer  melts.  Low  molecular  weight 

polymer  chains  move  essentially  independently  and  the  relaxation  times  of 

these  chains  scale  as  M^.  As  the  molecular  weight  increases,  entanglements 

with  neighboring  chains  are  more  frequent  and  large  scale  cooperative  motions 

couple  via  the  entanglement  junctions  to  the  CS  of  the  polymer  system.  The 

longest  relaxation  times  of  the  chain  then  scale,  by  Eq.  (4),  as  n^' 

If,  in  addition,  spatial  entanglements  are  imposed  by  confining  the  chain 

within  a  tube-like  region,  quantitatively  accurate  predictions  are  obtained 

concerning  the  viscosity  q  (where  qaM^),  the  recoverable  compliance  J°,  the 

plateau  modulus  G°  and  other  quantities  in  the  terminal  regime  for  both  linear 

and  branched  polymers.  For  linear  polymers,  n  is  typically  between  0.4  and 

0.5,  and  this  results  in  3.3<p<4  and  2<J°G°<3.  For  branched  polymers,  as 

determined  from  dispersion,  n  typically  has  larger  values  which  implies  larger 

M  and  J°G°.  These  quantities  can  simultaneously  be  correlated  with  changes  in 

the  temperature  dependence  by  the  flow  activation  energy.  All  these  features 

are  in  remarkable  agreement  with  data.  The  ability  of  this  model  to  describe 

13 

polymer  melt  dynamics  is  far  superior  to  the  reptation  model.  '  Details  of 

12 

the  polymer  melt  model  are  presented  elsewhere. 


Other  implications  of  Eq.  (4)  include  dielectric  and  structural  relaxa- 

14 

tions  near  for  polymers  and  small  molecule  glasses,  y  relaxation  of  bulk 
polycarbonate,*^  and  volume  recovery***  of  PVAC.  Thus  the  dispersion  given  by 
Eq.  (1)  and  the  shifted  relaxation  times  given  by  Eq.  (4)  play  important  roles 


in  a  large  number  of  materials  and  relaxation  phenomena. 


IX.  Calculations  Using  Time  Scale  Transformations 

Within  the  model  of  CS  excitations  of  Sec.  VIII,  the  time  scale  trans¬ 
formation  of  Sec.  II  becomes  0(t)=A  t*  n  where  A  =(l-n)  *e  n^ui  n.  The  use  of 

n  n  c 

this  technique  considerably  simplifies  calculations  of  relaxation  phenomena. 

This  is  illustrated  here  by  considering  hole  transport  across  a  biased  chal- 

cognide  glass  sample  of  thickness  L.  In  the  6-frame,  hopping  relaxation  at 

any  site  will  have  the  same  time  independent  transition  rate  Wg .  The  current 

in  the  6-frame,  1(0) 5  is  easily  determined  and  the  transient  current  in  the 

laboratory  t-frame  can  be  calculated  by  I(t)=I(0) (d8/dt) .  With  0(t)  given  as 

above,  it  is  found  that  I(t)at  n  for  t«t_,  and  a  more  rapid  drop-off  for 

t»t_.  The  transit  time  t„a(E/LW  )*^*  °  exhibits  superlinear  thickness  and 
i  -.is 

electric  field  dependences  and  is  thermally  activated  with  energy  E*=E^/(l-n) 

if  W  is  thermally  activated  with  energy  E. . * 

8  A 


X.  •  Applicability  of  Fractional  Exponential  Relaxation 

The  fractional  exponential  relaxation  function  given  by  Eq.  (1)  along . 
with  the  shifted  relaxation  times  of  Eq.  (4)  lead  to  unifying  predictions  in 
many  fields  and  materials.  These  include  dielectric  relaxation,  viscoelastic 
relaxation,  nuclear  spin  relaxation,  flicker  1/f  noise,  generation- 
recombination  noise,  transient  electrical  transport,  transient  capacitance, 
photoluminescence,  volume  and  enthalpy  recovery,  physical  aging,  plasticity 
and  yielding,  and  polymer  melt  dynamics.  Further  study  can  be  expected  both 


to  yield  deeper  insight  into  the  nature  of  low  frequency  fluctuation,  dissipa¬ 
tion  and  relaxation  phenomena  in  condensed  matter,  and  to  result  in  predic¬ 
tions  for  many  other  materials  and  applications. 


XI.  Summary  and  Conclusions 

The  dominance  of  fractional  time  exponential  relaxation  of  the  form 
Eq.  (1)  in  diverse  phenomena  has  been  emphasized  here.  Such  relaxation  func¬ 
tions  require  that  the  accompanying  transition  rates  have  power  law  time 
dependences.  These  dependences  can  be  physically  motivated  in  several  dif¬ 
ferent  ways.  However  the  dispersion  of  the  relaxation  predicted  b  q.  (1)  is 
observed  to  correlate  with  several  physical  quantities  including  -tivation 
energies,  and  temperature  and  molecular  weight  dependences.  A  pb  model 
in'  which  a  relaxing  species  interacts  with  a  set  of  low  energy  states  is  able 
to  predict  such  correlations.  It  is  important  to  stress  that  a  proper 
description  of  relaxation  requires  not  only  fractional  powers  of  time  in  the 
dispersion  as  in  Eq.  (1),  but  also  relaxation  times  which  shift  with  changes 
in  dispersion  as  in  Eq.  (A). 
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